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M  ABSTRACT  (MAXIMUM  200  WORDS) 


A  supersonic,  aerodynamic  computational  model,  which  is  the  basis  of  the  NANC  code, 
has  been  extended  to  compute  dynamic  derivatives.  The  extension  is  to  the  inviscid  conti  i- 
bution  of  constant  angular  rates  and  axial  accelerations. 

The  body  geometry  limitations  are  the  same  as  for  the  steady-state  model.  Here,  a 
pointed  body  or  equivalent  pointed  body  is  assumed  for  low  Mach  numbers;  at  higher  Mach 
numbers,  the  effect  of  axial  acceleration  is  neglected  and  the  body  may  be  blunt.  The  body 
may  be  noncircular  with  planar  discontinuities,  including  inlets,  with  fins  (up  to  six  per  fin 
set),  which  lie  on  a  cylindrical  coordinate  ray. 

For  the  low  Mach  number  range,  the  original  second-order  potential  model  has  been 
extended  for  angular  rate  derivative  prediction.  For  the  acceleration  rate  derivatives,  a 
"hybrid”  first-  and  second-order  model  has  been  developed. 

For  the  high  Mach  number  range,  an  equivalent  angle-of-attack  vector  is  defined  and 
combined  with  local  solution  models. 

Computational  comparisons  are  made  with  experimental  data,  primarily  for  pitch  and 
roll  damping  derivatives.  -  '  ■ 
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FOREWORD 


This  work  represents  an  extension  of  the  work  reported  in  Naval  Surface 
Warfare  Center  (NSWC)  TR  86-253.  The  latter  report  presented  computational 
methods  for  predicting  aerodynamic  loading  for  supersonic  Mach  numbers  on  non- 
axisymmetric  flight  vehicles  at  constant  or  steady  incidence.  The  extension  is  for  the 
prediction  of  aerodynamic  loading  or  dynamic  derivatives  associated  with  constant 
body  axis  angular  rates  and/or  acceleration.  The  resulting  computer  program  allows 
one  to  predict  roll  and  pitch  damping,  fin  Magnus,  and  other  dynamic  derivatives  for 
the  preliminary  and  intermediate  design  stage. 

Support  for  the  work  was  provided  by  the  following  sponsors: 

1.  Aerodynamics  and  Structures  Block  of  the  Surface-Launched  Weaponry  Tech¬ 
nology  Program  Project  Numbers  RA11G13/RU11G14. 

2.  Air-Launched  Antisurface  Weaponry  Technology  Project  Number  RU 11811. 

This  report  was  reviewed  and  approved  by  Dr.  T.  J.  Rice,  Head,  Aeromechanics 
Branch  and  C.  A.  Cooper,  Head,  Missile  Systems  Division. 
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1.0  INTRODUCTION 


Preliminary  design  requires  the  estimate  of  aerodynamics  for  a  large  set  of 
free-stream  and  geometric  variations.  To  keep  computational  costs  reasonable, 
rapid  (but  reasonably  accurate)  methods  are  sought  and  utilized. 

For  cruciform  finned,  axisymmetric  bodies  at  low  angle  of  attack,  component 
buildup  methods1.2  can  be  utilized  for  a  restricted  range  of  configurations  and  free- 
stream  conditions  that  are  applicable  to  current  designs. 

Linear  surface  singularity  methods3-4  have  been  highly  developed  for  complex 
configurations,  including  high  angle-of-attack  vortex  modelling.  Set-up  and  run 
times  are  fairly  long. 

The  second-order  Van  Dyke  model5,  which  corrects  for  compressibility  effects, 
was  extensively  modified  and  adapted  for  noncircular  bodies  with  planar  discon¬ 
tinuities,  including  inlets,  and  multi  sets  of  fins.6  9  The  NANC  code1912  includes 
the  second-order  potential  model  plus  a  local  solution  model  based  on  the  methods  of 
Reference  13. 

Dynamic  derivative  estimates  are  calculated  in  the  NSWC  Aeroprediction  code 
of  Reference  1.  Body-alone  contributions  are  given  by  an  empirical  data  fit.  Fin- 
alone  contributions  are  given  by  thin-wing  theory.  Interference  modelling  is  incom¬ 
plete  when  compared  with  the  static  case.  Thin-wing  theory  is  based  on  conical 
solutions  that  are  not  valid  for  small  aspect  ratio  fins. 

Dynamic  derivative  estimation  is  primarily  of  interest  for  unguided  appli¬ 
cations.  Pitch  damping  is  usually  ignored  for  guided  applications  since  the  autopilot 
provides  control  surface  moments,  which  are  proportional  to  the  pitch  rate  and 
significantly  larger. 

Classical  unsteady  aerodynamics  is  based  upon  the  unsteady  linear  potential 
equation.14  The  unsteady  potential  equation  may  be  applied  to  harmonic  analysis  of 
the  rigid  motion  or  to  aeroelastic  applications.  The  harmonic  gradient  method  of 
Reference  15  is  a  singularity  collocation  code  utilizing  a  complex  potential.  The  total 
pitch  damping  for  low  frequency  is  the  same  as  for  the  constant  pitch  and  accel¬ 
eration  rate  estimate. 

The  NANC  code  is  modified  in  this  report  to  compute  forces  and  moments 
which  are  a  function  of  constant  axial  rotation  rates  and  moment  center  accel¬ 
erations.  The  remainder  of  the  report  assumes  that  the  reader  is  familiar  with  the 
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earlier  work  summarized  in  Reference  10  and  concentrates  on  the  dynamic 
derivative  modifications. 


2.0  GEOMETRY,  KREE-STREAM  VELOCITY 
AND  FORCE  CONVENTIONS 


Assumptions  concerning  body  and  fin  geometry  are  the  same  as  in  Reference 
10.  Figures  1  through  3  are  taken  from  Reference  10.  However,  the  body  is  assumed 
to  have  a  pointed  nose  for  the  low  Mach  number  range. 

A  static  blunt  body  model  establishes  a  matching  plane  of  velocity  component 
data  by  matching  a  modified  Newtonian  pressure  distribution  and  utilizing  conical 
potential  functions  and  other  assumptions.  A  blunt  body  model  for  the  plunging 
acceleration  problem  was  not  deemed  feasible.  A  model  for  pure  axial  rotation  is 
only  somewhat  more  feasible. 

Computations  utilize  a  cylindrical  coordinate  system  and  the  "thin-fin”  approx¬ 
imation  with  the  fin  midplane  lying  on  a  cylindrical  coordinate  ray.  The  body  is 
divided  into  sections  by  planar  discontinuities  (including  planar  inlets).  Any  section 
except  the  first  may  have  a  set  of  fins  with  up  to  six  fins.  See  Reference  10  for  a  more 
detailed  discussion. 

Force  convention  is  the  same  as  used  in  earlier  work.  The  axial  force,  Fa,  acts 
in  the  x  direction;  the  normal  force,  F.\,  acts  in  the  z  direction;  and  the  side  force,  Fy, 
acts  in  the  negative  y  direction.  The  roll  moment,  Mf,  acts  in  the  negative  x 
direction;  the  yawing  moment,  Mn,  acts  in  the  negative  z  direction;  and  the  pitching 
moment,  Mm,  acts  in  the  negative  y  direction. 

The  main  difference  between  the  current  work  and  earlier  work  is  the  body- 
axis  oriented  equivalent  dimensionless  free-stream  velocity  vector.  In  Cartesian 
coordinates,  the  velocity  vector  relative  to  body  axes,  in  dimensionless  form,  is  given 
as 


V  -q'z  +  r'y 

—  +  -  -I-  cos«  cosp 

V  V 


+ 


—  r'(x  —  x')  +  p'z 
V 


+  sinP'  j' 


(2-1) 


+ 


p'y  4-  q'(x  —  x') 
V 


+  sin«  cosP’ 


k' 


<L  f  <1,' 


lim  t  — »  0 
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Note  that  the  time  derivatives  of  the  axial  rotation  rates  could  be  included  as  well, 
but  are  neglected.  p\  q',  r'  are  the  axial  rotation  rates  in  the  negative  x,  y,  and  z 
directions,  respectively,  and  uo,  vo,  and  wo  are  the  corresponding  acceleration  rates. 
For  the  total  pitch  damping  problem,  q’  and  wo  are  the  pertinent  input  terms.  «  and 
P'  are  the  angle  of  attack  and  side  slip,  respectively. 

An  additional  symmetry  operation  variation  for  LM  (defined  in  Reference  10)  is 
introduced  here.  LM  =  1  is  for  axial  force  as  the  only  nonzero  force  term.  LM  =  2  is 
for  axial  force  and  rolling  moment  as  the  only  nonzero  forces.  The  new  symmetry 
value  for  LM  =  2  is  for  a  circular  body  with  quarter-plane  symmetric  fin  distribution 
and  a  roll  control  deflection  or  constant  roll  rate.  For  the  LM  =  2  case,  the 
computation  is  LM  =  1  before  the  first  fin  set  and  LM  =  2  subsequently.  LM  =  3  is 
for  the  axial  force,  normal  force,  and  pitching  moment  nonzero.  LM  =  5  is  for  axial 
force,  side  force,  and  side  moment  nonzero.  LM  =  6  is  for  a  full  plane  computation. 
LM  is  defined  for  each  section  of  the  body  and  must  increase  or  remain  the  same  as 
the  computation  is  marched  down  the  body. 


3.0  THEORETICAL  DEVELOPMENT 


3.1  FIRST-  AND  SECOND-ORDER  POTENTIAL  EQUATIONS 
The  full  nonlinear,  dimensionless,  potential  equation  is 


a~v  •  Q  --  m2 


j_  f4>  J_  Q  VQ2 
v7  at2  +  v,  at  2 

r 


(3-1) 


Here,  a  is  the  speed  of  sound  divided  by  and  Q  =  V(J>  is  the  dimensionless  velocity 
vector.  Equation  (3-1)  is  appropriate  to  a  body  moving  through  a  fluid  at  rest.  It  is 
derived  from  the  continuity  equation,  momentum  equations,  and  the  Bernoulli 
relationship.  For  body-fixed  coordinates, 


a 

at 


a 

at 


v-v 


(3-2) 


V  = 


-V 


(3-3) 


In  addition,  the  relative  velocity,  qr  =  Q  -  V,  is  introduced.  The  continuity 
equation  then  becomes 


l  l  ap  i  /  o  a  \ 

- —  +  V  •  <|  +  -  (q  -V)p  =  0  (3-4) 

vr  p  at  ’  P  r 


3 


NSWC  TR  89-99 


The  Bernoulli  relationship  for  body  coordinates  becomes 


M- 


o  [  1  # 


q;  -  v- 


v  at 


(3-5) 


V  -  1  \  -  1 


The  density  may  be  eliminated  from  Equation  (3-4)  by  using  Equation  (3-5)  to  yield 


a2V  -q  =  \r 
1  r  * 


•  V 


fi,  “  v  >  l  a‘4>  a  /  v  v  aQ 

-  + - —  +  —  -  +  q  •  — 

2  1  v2  at2  *  '  2  '  <n 


dQ  1 

at  f 


(3-6) 


a“  -  1  + 


V  -  1 


M- 


V“  -  q" 


cxb 

2<  — ) 
at 


(3-7) 


A  small  disturbance  from  the  static  free-stream  dimensionless  velocity  vector  is 
given  by 


qt.  =  f  +  <q(>r  -  «  )  +  q  =  i'  +  q>  q  (3'8) 

Substitution  of  Equation  (3-8)  into  Equation  (3-6)  and  neglecting  higher  order 
terms  yields  the  first-order  wave  equation 


V-q', 


r 


(3-9) 


q't  is  the  first-order  velocity  vector,  u'j  the  x  component  of  qj,  and  4>j  the  corresponding 
potential  function. 

The  first-order  problem  may  be  further  broken  down  into  three  steady  prob¬ 


lems. 


M2 


d>'  =  d>,  +  —  d>  +  V 

M  M  p2  Ta  3 


t  - 


xM 

'» 

p2v 


(3-10) 


Here  4>i  is  the  potential  function  associated  with  the  angle-of-attack  problem  and 
axial  rotation  rates.  <J>a  and  $b  are  equivalent  steady  potential  functions  associated 
with  the  plunging  rate  problem.  Equation  (3-10)  satisfies  Equation  (3-9)  when  <J>i, 
<j)a,  and  4>b  satisfy  the  steady  first-order  wave  equations 


(V 


Vo 


c?X 


(3-11) 
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(V  x  q .)  =  0  (3-12) 

f  c 

£  —  1,  a,  b;  c  stands  for  the  crossflow  plane  and  components. 

As  in  earlier  work,  an  improved  solution  to  the  first-order  pn  blem  is  obtained 
by  evaluating  the  neglected  nonlinear  terms  using  the  first-order  solution  and 
solving  a  nonhomogeneous  potential  equation.  Only  time-independent  terms  are 
considered. 


<3w 

(V  •  q  ,)  -  B2  — 1  =  M2  V  •  12 
*  dx 

(3-13) 

(V  x  q  )  =0 

it 

(3-14) 

12  =  12  i'  +  12 

\  C 

(3-15) 

n  1  >  (2  -  y)  4  2  (2  -  V )  2  ■>  , 

Qx~ui  2qi“  6  M.u.  -  2 

Y  —  1  , •>  2 

+  2  M;u;  +  q,  +  <vqt 

(3-16) 

•> 

qi  2  -  v  •>  •>  9 

Q,  -q.c  2  +  2  M;U>  +q*v<q.  +q^’ql) 

M2  q2 

I  '  to 

+  —  (v  —  1 )  —  q  +  q ,  X  (q  x  q  ) 

(5 

(3-17) 

Boundary  conditions  for  first-  and  second-order  problems  are 

<q  +  q  )  ■  n  =  0  (3-18) 

£  =  1,2 

The  equivalent  steady  problem  boundary  conditions  for  tiie  q^  problems  are 
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Here,  t)r,  /dx  and  tx  are  body  and  fin  slopes,  respectively. 


3.2  COMPUTATIONAL  COORDINATES  AND  GRIDS 

The  region  between  the  body  and  the  Mach  cone  is  mapped  to  the  rectangular 
region  shown  in  Figure  2  by  the  shearing  transformation 


Further  clustering  transformations  E,  =  HO  and  <t>  =  are  used  for  nonuniform 
gridding  of  the  B  and  £,  variables.  The  functional  dependencies  are  not  given 
explicitly.  See  Section  3.2  in  Reference  10  for  more  details. 


3.3  NUMERICAL  METHODS 

Most  of  the  numerical  methods  are  as  reported  in  Reference  10.  However,  some 
of  them  have  been  changed. 

The  implicit  formulation  for  the  first  body  section  is  as  in  earlier  work.  The 
velocity  vector  advancement  Equations  (3-11)  through  (3-14)  are  single  second-order 
equations  in  a  potential,  $  =  xF.f."  F  is  known  as  a  conical  potential  function. 
However,  actual  conical  similarity  for  the  total  pitch  damping  problem  requires  the 
functional  form  of  <£  =  x  (G  +  xH).  The  pressure  distribution  on  a  cone  is  now  linear 
with  x  instead  of  constant.  Therefore,  the  computation  uses  more  than  two  steps  for  a 
solution  for  a  cone.  The  solution  at  x  =  0  is  assumed  to  be  conical.  However,  the 
resultant  computation  for  a  numerical  solution  using  a  few  marching  steps  does  not 
vary  significantly  from  a  solution  developed  based  on  the  true  conical  similarity  of  <J) 
=  x  (G  +  xH)  and  two  marching  steps. 

At  body  planar  discontinuities  and  supersonic  leading  and  trailing  edges,  the 
jump  in  various  velocity  components  is  obtained  by  application  of  the  method  of 
"weak  solutions”1  r’  combined  with  a  downstream  boundary  condition  and  conser¬ 
vation  of  the  velocity  component  tangent  to  the  discontinuity  edge  condition.  At  a 
subsonic  leading  edge,  the  solution  for  the  velocity  jumps  does  not  exist  and  the 
conservation  relationship  provided  by  the  method  of  "weak  solutions”  must  be 
replaced  by  a  heuristic  one.  The  "weak  solution”  conservation  relation  is  modified 
due  to  the  modified  Equation  (3-13). 

For  all  sections,  except  the  first,  the  conservation  velocity  vector  advancement 
equations  are  solved  using  a  MacCormack'7  predictor-corrector  scheme  for  points  not 
on  a  solid  surface.  Body  and  fin  surface  velocities  are  advanced  by  characteristic 
compatibility  relations  (two  in  number)  combined  with  a  solid  surface  boundary 


6 


NSWC  TR  89-99 


condition. 18  At  a  corner  where  the  fin  and  body  meet,  no  unique  solution  is  possible. 
The  advancement  equations  for  the  fin-body  junction  have  been  modified.  The  fin- 
body  junction  line  is  assumed  to  be  a  streamline.  Therefore,  fin  and  body  boundary 
conditions  and  one  other  condition  are  needed  to  solve  for  the  three  velocity 
components. 

One  estimate  of  the  axial  velocity  component,  u,  can  be  obtained  from  one  body 
advancement  quantity  and  the  body  boundary  condition. 


=  ubP  v/)  +e0  + vb  -£0wb 

(3-22) 

fJY|, 

-cowb=  —  (uB  +  uj-v»+£flw; 

(3-23) 

Here,  u,  v,  and  w  are  the  axial,  radial,  and  0  components  of  velocity,  ey  =  l/rb(drb/d0). 
The  B  subscript  stands  for  the  body.  For  the  first-  and  second-order  problem,  ui,  v;, 
and  w'rare  free-stream  velocity  components.  For  the  <J>b  problem,  q’x=  q*.  For  the  (f>a 
problem,  u'r  =  urx/V„  -  4>b,  v;  =  v*x/Vr,  and  w^  =  Wxx/V*.  uB  is  obtained  from  Equa¬ 
tions  (3-22)  and  (3-23).  Eb  is  known. 

A  second  estimate  of  the  axial  velocity  component  is  obtained  from  the  fin 
advancement  quantity  and  the  fin  boundary  condition 


Ef  =  ±  PuF  +  wf. 


WF  =  -  wj  +  MuF  +  u'r\ 


(3-24) 

(3-25) 


Here,  Ep  is  known,  up  may  be  obtained  from  Equations  (3-24)  and  (3-25).  The  upper 
sign  is  for  the  0  >  0fm  side  of  the  fin. 

The  final  values  of  u,  v,  w  for  the  fin-body  junction  are  obtained  from 

u  =  \  (UB  +  V  (3’26) 


v  -  cw  =  —  (u  +  u’„)  -  v;+  c  w; 

0 


w  =  -  w'  +  t  (u  +  u' ) 

r  X 


(3-27) 

(3-28) 


Note  that  $b  need  only  be  determined  on  the  solid  surfaces  and  is  needed  for  the  <t>a 
problem  boundary  condition  and  evaluation  of  the  Bernoulli  pressure  coefficient 
relationship.  An  advancement  equation  for  <J>b  is  given  by 


—  <  1  -  £,)  +  -  v  +  t  w 

rlx  (1  b  x  h 


(3-29) 
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On  a  body  surface  or  interior  point,  tx  is  set  to  zero. 

A  subsonic  leading  edge  requires  a  modification  of  the  jump  relations.  For  all 
cases, 


k.j 


k-l.j 


2  2 
r0  ~  rkj 


(3-30) 


k  is  the  grid  index  in  the  r  direction;  j  is  the  index  in  the  0  direction.  Equation  (3-30) 
combined  with  downstream  boundary  condition  and  conservation  of  the  velocity 
component  tangent  to  the  edge  provides  a  solution  for  the  downstream  values  of 
velocity.  The  square  root  ratio  is  limited  to  a  value  of  2.  ro  is  the  radius  of  the 
leading  edge.  Equation  (3-30)  has  the  well  known  square  root  singularity  for  a 
subsonic  leading  edge. 


3.4  LOADING  COEFFICIENTS  AND  OTHER  NUMERICAL  CONSIDERATIONS 

The  axial  acceleration  rates  contribution  to  velocities  in  the  limit  as  t  -*  0  is 
given  by  differentiation  of  the  last  two  terms  of  Equation  (3-10). 


q. 


ix 


P2 


M2 

t 

P2 


(3-31) 


The  first-  or  second-order  pressure  coefficient  is  then  given  by 


Cp  =!j3-5  -  1|/(.7M;) 

(3-32) 

jf=l+.2M>2-Q2f-2<V 

(3-33) 

=  V  +  (U  +  q<’ 

(3-34) 

e  =  1,2 

This  is  a  "hybrid”  model  for  the  plunging  rate  case  since  qax  is  a  first-order  potential 
quantity. 

The  inviscid  loading  coefficients  are  as  given  in  Section  5.0  of  Reference  10. 

Smoothing  of  the  MacCormack  vector  quantities  is  an  input  option  for  sections 
with  fins.  It  is  particularly  needed  for  subsonic  leading  edges.  Q|,  with  qax  =  0  is 
used  as  a  weighting  function  for  first-  and  second-order  vector  terms  as  in  Section 
3.11  of  Reference  10.  Q|,  with  qux  *  0  is  used  as  a  weighting  function  for  the  $a  and 
4>b  problem  vector  terms. 
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The  explicit  marching  solution  will  fail  when  X  becomes  negative  in  Equation 
(3-32).  For  q,B  =  0,  X  is  set  to  zero  or  the  vacuum  value  and  the  velocities  are  adjusted 
as  discussed  in  Section  3.11  of  Reference  10.  For  q*  *  0  and  X  negative,  the  solution 
is  halted  and  an  error  message  is  written  to  the  output  File. 


3.5  HIGH  MACH  NUMBER  SOLUTION 

The  potential  model  breaks  down  at  Mach  numbers  where  the  origin  Mach  cone 
crosses  the  body  surface,  q,.  is  set  to  zero  for  this  local  solution  model. 

The  velocity  vector,  q*,  is  normalized  to  1  as 

(3-35) 

q^  is  then  used  to  find  the  turning  angle  between  an  effective  free-stream  velocity 
vector  and  the  local  surface  normal.  These  are  combined  with  the  local  solution 
methods  of  Reference  13.  It  is  expected  that  these  local  solution  methods  will  be 
inaccurate  below  M,.  =  4.  Note  that  blunt  bodies  may  be  considered. 


4.0  EVALUATION  OF  THE  NUMERICAL  METHODS 


Evaluation  is  almost  exclusively  made  by  comparison  with  experimental  data. 
The  data  sources  are  ballistic  range  data  or  wind  tunnel  tests.  Dynamic  derivatives 
must  be  based  on  parameter  estimation  techniques  associated  with  kinematic  data 
and  an  assumed  aerodynamic  model.  In  general,  the  accuracy  and  repeatibility  of 
the  estimated  coefficients  is  much  worse  than  for  static  coefficients. 


The  code  development  is  for  nonaxisymmetric  bodies  with  combined  angle  of 
attack,  control  deflection,  rotation  rates,  and  axial  accelerations.  Data,  however,  is 
mostly  available  for  Cipr  and  Cmq'  +  Cmi  for  circular  bodies  at  zero  incidence. 
Computations  are  separate  for  Cip',  Cmq',  and  Cmd. 

The  dynamic  derivatives  are  given  in  general  form  as 


C  = 


M 


^KXR<^D(a>iXD/  V®) 


(4-1) 


Sr  is  a  reference  area;  xr  is  a  reference  length;  M;  is  a  force  or  moment;  Qd  is  the 
dynamic  pressure;  xo  is  a  length  associated  with  a  reduced  frequency;  and  co;  is  an 
axial  roll  rate  or  dimensionless  acceleration. 
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4.1  BODY- ALONE  COMPARISONS 

The  reference  area  for  bodies  will  be  the  maximum  cross-sectional  area. 

The  first  computation  is  for  a  2.98  caliber  length  cone.  The  moment  center  is 
located  2.18  calibers  from  the  nose,  xr  is  the  maximum  body  diameter  and  x[)  is  the 
maximum  body  radius.  Unless  otherwise  indicated,  references  will  not  vary  in  the 
rest  of  the  report.  Figure  4  compares  Cmq'  +  Cm*  computational  values  with  values 
extracted  from  the  ballistic  range  data  of  Reference  19.  The  usual  lack  of  repeat- 
ibility  of  the  data  is  shown  due  to  different  initial  yaws  and,  hence,  epicyclic  history. 
Figure  5  compares  Cmq'  +  Cmi  computational  values  for  a  cone-cylinder  and  the 
ballistic  range  data  of  Reference  19. 

The  next  comparison  computations  are  for  the  Army-Navy  Spinner  config¬ 
urations.  These  configurations  have  2-caliber  secant  ogive  noses  with  arc  radii  twice 
that  of  a  tangent  ogive  and  various  body  lengths  and  moment  center  locations. 
Figure  6  compares  Cmq’  +  Cmi  computations  with  data  for  a  5-caliber  length  body. 
Note  the  great  differences  between  different  range  tests  and  wind  tunnel  tests.  Also 
plotted  are  Cm«  and  a  slender  body  value.  Figure  7  compares  Cmq'  +  Cmi  compu¬ 
tations  with  data  for  a  9-caliber  length  body.  GE-Spinner  refers  to  the  empirical 
curve  fit  of  Reference  20. 

The  final  body  computational  example  is  taken  from  Reference  15.  xr  and  xd 
are  cone  lengths,  L,  for  an  elliptic  cone.  Semi-minor  to  semi-major  axis  ratio  is  a/b  = 
.75.  a/L  =  .0866,  x’  =  0.  Figure  8  compares  C\q'  +  C\i  first-  and  second-order 
computations  with  that  of  Reference  15.  Figure  9  compares  Cmq'  +  Cm£.  Compu¬ 
tations  for  Cvi  show  significant  differences  between  first-  and  second-order  order. 
The  computation  of  Reference  15  is  a  first-order  computation  using  different  numer¬ 
ical  methods.  Three  computational  planes  seems  to  be  adequate  for  this  case. 


4.2  BODY-TAIL  CONFIGURATIONS 

The  first  configuration  considered  for  this  section  is  the  Basic  Finner  of  Figure 
10.21 ,22  Figure  1 1  shows  a  comparison  of  data  with  computation  for  the  roll  damping 
derivative.  The  loading  for  the  roll  damping  problem  is  low  near  the  body  and 
increases  with  span  distance.  The  roll  moment  loading  increases  even  more  rapidly 
with  span  distance.  At  a  side  edge,  the  loading  drops  to  zero.  Load  integration 
routines  were  modified  to  try  to  account  for  the  sharp  drop  off  of  loading  near  a  side 
edge.  This  improved  C|p'  prediction  at  lower  Mach  numbers.  Figure  12  depicts  a 
Magnus  moment  derivative,  C.\p -r,  computation  comparison  with  data.  A  body- 
alone  estimate  using  the  empirical  computation  of  Reference  20  indicates  a  variation 
of -.4  to  3  for  the  Mach  number  range  of  1.2  to  3.  The  computation  here  only  follows 
for  the  data  trend  but  not  the  magnitude.  Figure  13  depicts  a  total  pitch  damping 
comparison.  A  trend  of  being  above  the  data  at  lower  Mach  numbers  is  indicated 
here. 
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The  next  configuration  is  that  of  Figure  14(b)  taken  from  Reference  23.  The 
equation  of  the  body  in  calibers  is 

rb=  -511  -  (1  -  x/6.25)2!3'4  0  <  x  <  6.25  (4-2) 

=  .5,  x  >  6.25 

Since  drb/dx  and  d2rb/dx2  are  singular  at  x  =  0,  the  nose  is  approximated  as 


r  =  .297 1 6844x  -  .13 167349 lx"  +  .034340881 x3  0  <  X  <  1  (4-3) 

b 


Equation  (4-3)  matches  zero  through  second  derivatives  of  Equation  (4-2)  at  x  =  1. 
The  total  pitch  damping  comparison  is  shown  in  Figure  15.  Reference  area  is  the 
extended  to  centerline  fin  planform  area  (two  fins).  Reference  xr  =  2xj)  =  the  mean 
aerodynamic  chord, 


„  <C2  +  C2  +  C  C  ) 

2  r  t  t  r 

3  C  +  C 

t  r 


(4-4) 


Here,  Cr  and  Ct  are  root  and  tip  chords  of  the  extended  fin.  The  moment  center  is  .2C' 
from  the  apex  of  the  wing  extended  to  the  body  centerline.  Leading  and  trailing 
edges  are  supersonic  and  the  computation  compares  fairly  well  with  the  data. 
Computations  were  adequate  for  the  configuration  of  Figure  14(a)  for  supersonic 
leading  edges,  but  very  poor  for  subsonic  edges  when  the  moment  center  is  .35C'  from 
the  extended  wing  apex.  The  reason  is  the  inaccurate  computation  of  Cmi.  For 
moment  centers  which  are,  more  typically,  not  close  to  the  point  of  action  of  Cx«,  the 
sensitivity  and  accuracy  is  better. 

The  next  configuration  considered  is  depicted  in  Figure  16.  A  Cmq’  +  C  mi 
comparison  for  this  configuration  is  shown  in  Figure  17.  Here,  the  trend  of  over¬ 
prediction  is  reversed. 

The  next  body-tail  configuration  considered  is  the  fiechette  of  Figure  18  taken 
from  Reference  24.  Reference  24  (unpublished)  compares  one  of  the  routines  for 
computing  total  pitch  damping  in  Reference  1  with  available  data  for  a  large  number 
of  body-alone  and  body-tail  configurations.  The  XM-144  body  end  radius  of  .005 
calibers  is  approximated  as  .05  calibers.  Figure  19  compares  computations  with  data. 
For  this  case,  the  leading  edge  of  the  fins  is  subsonic  below  Mach  numbers  of  3.25. 
The  data  was  surprisingly  smooth  for  this  configuration.  The  boat  tail  has  a 
significant  effect  on  the  computational  results. 

The  final  body-tail  configurations  0f  Figure  20  has  six  fins  and  a  boat-tail 
angle  of  2  deg.  Figure  21  compares  C|p'  potential  and  local  computations  with  the 
PNS  computations  of  Reference  25.  The  sharp  rise  close  to  sonic  leading  edge 
conditions  is  predicted  by  the  thin-fin,  fin-alone  methods  of  Reference  1.  The  boat  tail 
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significantly  affects  the  solution.  Figure  22  shows  a  pitch  damping  computational 
comparison  with  data  for  a  configuration  close  to  that  of  Figure  20.  The  configuration 
and  data  are  from  Reference  24.  Here,  the  conical  nose  length  is  3.32  calibers  and  the 
moment  center  is  7.14  calibers  from  the  nose.  From  Figures  21  and  22,  one  can  see 
that  the  local  solution  is  of  limited  value. 


4.3  BODY-WING-TAIL  OR  BODY-CANARD-TAIL  CONFIGURATIONS 

Here,  the  additional  interference  of  forward  lifting  surfaces  on  the  tail  is  the 
primary  phenomenon.  As  in  earlier  work,  it  is  assumed  that  the  velocity  vector 
downstream  of  a  supersonic  trailing  edge  lies  on  a  constant  cylindrical  ray  plane,  0 

=  0f. 


The  first  computational  comparison  is  for  the  B-C-T  configuration  of  Figure  23. 
Figure  24  shows  a  computational  comparison  with  experimental  data.  The  trend  at 
lower  Mach  numbers  is  as  noted  earlier.  Here,  the  carryover  canard  to  tail  of  Cmi 
grows  as  the  Mach  number  decreases.  Evaluation  of  the  plunging  acceleration  con¬ 
tribution  of  Equation  (3-31)  involves  Ml/02  and  differences  which  probably  become 
sensitive  as  the  Mach  number  is  decreased.  Computations  of  Cip’  aiso  indicate  a 
carryover  effect  which  is  not  accounted  for  by  component  superposition  methods.  At 
Mo,  =  1.76,  C|p'  =  -199.3  for  the  total  configuration.  For  the  same  Mach  number,  the 
total  of  individual  canard  and  tail  configurations  adds  up  to  Cip  =  -171.1. 

Figure  25  compares  total  pitch  damping  with  data  and  a  semi-empirical  esti¬ 
mate  for  the  RFL  122  configuration  of  Reference  26.  The  ordinate  scale  was  not  given 
and  is  inferred  here.  The  Cma  carryover  phenomenon  seems  to  be  less  extreme  in  this 
case.  The  semi-empirical  estimate  seems  to  be  worse,  assuming  the  data  fit  is  correct. 
A  Cip’  computation  at  M<r,  =  1.5  yields  -89.5  for  the  total  configuration.  The  wing 
contribution  to  C|p  at  the  same  Mach  number  is  -94.0. 

The  final  computation  is  for  the  B-W-T  configuration  of  Reference  27.  Total 
pitch  damping  comparison  and  configuration  are  depicted  in  Figure  26.  Geometry 
and  the  moment  center  location  are  given  in  Reference  27.  The  geometry  was  scaled 
from  Reference  27.  x'  is  assumed  to  be  at  x'  =  8.5  calibers  from  the  nose.  The  Cmi 
carryover  phenomenon  is  noted  as  before.  Missile  refers  to  an  Aerospatiale  semi- 
empirical  code. 
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5.0 CONCLUDING  REMARKS 


The  computational  methods,  developed  for  computing  static  aerodynamic 
coefficients  for  noncircular  bodies  at  supersonic  Mach  numbers,  have  been  extended 
to  the  computation  of  dynamic  derivatives. 

Geometric  limitations  are  the  same  as  for  the  earlier  work.  However,  it  is 
assumed  that  the  bodies  are  pointed. 

The  weakest  elements  in  the  original  code  are  exaggerated  here  since  moments 
are  almost  exclusively  computed.  The  variation  of  pressures  close  to  subsonic  lead¬ 
ing  and  side  edges  is  affected  by  edge  singularities.  C|p'  prediction  is  poorest  for  fins 
with  low  aspect  ratios.  Cmq'  +  Cmi  prediction  is  poorest  when  the  leading  edge  is 
subsonic  and  the  point  of  action  of  C\(,'  +  C\i  is  close  to  the  moment  center.  The 
carryover  effect  for  Cm«  is  overpredicted  when  the  leading  edges  are  subsonic  for 
bodies  with  two  sets  of  lifting  surfaces. 

The  methods  developed  are  capable  of  computing  dynamic  derivatives  that  are 
usually  not  computed  or  measured  for  in-plane  and  out-of-plane  cases. 

Computational  comparison  with  data  is  fairly  good  for  most  cases.  In  general, 
the  repeatibility  of  dynamic  derivative  experimental  values  is  much  poorer  than  for 
static  derivatives. 

The  current  computer  code,  implementing  the  dynamic  derivative  models,  re¬ 
quires  about  330,000  octal  storage  locations  for  a  15-by-60  grid  for  full  plane 
computations.  Computational  times  for  axial  rate  computations  are  about  the  same 
as  for  static  computations.  Plunging  rate  computations  are  about  50  percent  longer. 
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FIGURE  10.  BASIC  FINNER  CONFIGURATION 


20 


-  CNPA  .  CLP 


NSWC  TR  89-99 


FIGURE  11.  BASIC  FINNERCip’ COMPARISON 


FIGURE  12  CMp  MAGNUS  DERIVATIVE  FOR  THE  BASIC  FINNER 


21 


NSWC  TE  89-99 


FIGURE  13.  TOTAL  PITCH  DAMPING  COMPARISON 
FOR  THE  BASIC  FINNER 


(a)  A  ■  3  SWEPT  WING  (3%  BICONVEX  CIRCULAR-ARC  SECTIONS)  A  “  3  UNSWEPT  WING  (3%  BICONVEX  CIRCULAR-ARC  SECTIONS) 


FIGURE  14.  ASPECT  RATIO  =  3  CONFIGURATIONS 
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FIGURE  21.  ROLL  DAMPING  COMPARISON  COMPUTATIONS 
FOR  THE  M735  PROJECTILE 
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FIGURE  22  TOTAL  PITCH  DAMPING  COMPARISON 
FOR  THE  M735  PROJECTILE 
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FIGURE  24.  TOTAL  PITCH  DAMPING 


FIGURE  25  TOTAL  PITCH  DAMPING  COMPARISON  FOR  THE  RFL  122 
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NOMENCLATURE 
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NOMENCLATURE 

a  Dimensionless  speed  of  sound 

aB  Reference  speed  of  sound  =  free-stream  value 

Ci  General  force  or  moment  coefficient 

Cp  Pressure  coefficient 

i  Grid  index  for  a  constant  x  plane 

i'  Unit  vector  in  x  direction 

j  Grid  index  for  constant  9  plane 

j'  Unit  vector  in  y  or  r  direction 

k  Grid  index  for  constants  plane 

k'  Unit  vector  in  z  or  0  direction 

LM  Aerodynamic  symmetry  mode 

Mj  Moment  or  force 

MB  Free-stream  Mach  number 

n  Unit  normal  vector  from  a  solid  surface 

p'  Rotational  rate  about  the  x  axis 

(NOTE:  All  q  velocities  are  nondimensional) 
q  Perturbation  velocity  vector 

q'  Rotation  rate  about  y  axis 

qa  Velocity  vector  associated  with  <{)a  potential 

qax  Velocity  vector  at  t  =  0  due  to  plunging  acceleration 
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9b 


9o» 


9r 


q,, 

q» 

q» 

9i 

«i 

q2 

9ic 

Q 

Qn 

QT? 

r 

r' 

rb 

r0 

Sr 

t 

tx 

u 

U„ 


Velocity  vector  associated  with  <J>b  potential 

Equivalent  free-stream  velocity  vector  relative  to  body  axes 

Q  -  V,  relative  velocity 

Equivalent  free-stream  velocity  vector  relative  to  body  coordinates  at  t 
=  0 

Free-stream  acceleration  relative  to  body  axes  (dimensionless) 

Equals  q0QO  -  i'  ^  " 

First-order  velocity  vector,  steady  part 
Total  first-order  velocity  vector 
Second-order  velocity  vector 
First-order  crossflow  vector 

Total  velocity  vector  of  fluid  in  body  axis  coordinates  relative  to  gas  at 
rest 

1/2  pB  V2,  dynamic  pressure 

Total  first-  or  second-order  velocity  vector,  t  =  0 

Radial  coordinate 

Rotation  rate  for  z  direction 

Body  radius 

Body  radius  to  outer  fin  edge 

Reference  area 

Time 

Fin  slope  in  x  direction 
x  component  of  perturbation  velocity 
x  component  of  c.  g.  acceleration 
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x  component  of 

r  component  of  perturbation  velocity 
y  component  of  c.  g.  acceleration 

Equals  -q0r,,  velocity  relative  to  gas  at  rest  in  body  axis  coordinates 
Magnitude  of  velocity  vector  of  c.  g.  at  t  =  0 
9  component  of  perturbation  velocity 
z  component  of  c.  g.  acceleration 

Cartesian  coordinate  from  nose  —  zero  incidence  direction 
Moment  center  distance  from  nose 
Reference  length 

Reference  length  associated  with  reduced  frequency 

Cartesian  coordinate  perpendicular  to  vertical  body  symmetry  plane  (to 
the  right  looking  downstream) 

Cartesian  coordinate  perpendicular  to  x,  y  plane  (up,  looking  down¬ 
stream) 

Angle  of  attack 

Equals  VmJ-  1 

Angle  of  side  slip 

Specific  heat  ratio 

Equals  (fJr^/fjB)/  r^ 

Equals  ^  (€,),  clustering  transformation  in  £,  direction 
Cylindrical  coordinate  angle  from  leeside  plane 
Equals  (r  -  r^/lx/p  -  ri,),  shearing  transformation 
Density 
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pr  Free-stream  density 

4>  Disturbance  potential 

(J)t  First-order  disturbance  potential  for  q,  =  0 

4>;  Equals  +  (Moo/p2)^.,  +  V^lt  -  xM2  /  (P2V00)]<|>b 

(Total  first-order  disturbance  potential  $a  and  4>b  are  equivalent  static 
problem  potentials.) 

(J>2  Second-order  disturbance  potential 

d>  Equals  (0)  -  clustering  transformation  in  B  direction 

X  Absolute  temperature  ratio 

ojj  General  term  in  computing  reduced  frequency 
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